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In the ELF range, the homogeneous, isotropic 
model ionosphere is not successful in explaining ob- 
served signal characteristics. This has prompted 
the introduction of gradually tapered ionosphere 
models by Wait [1960a, b, 1962] and Galejs [1961a, b, 
1962], which still have not considered the effects of 
ion motion and the earth's magnetic field. 

The purpose of this paper is primarily to esti- 
mate the probable effects on ELF propagation 
caused by the earth's magnetic field in combination 
with electron and ion densities which vary gradu- 
ally with height. To this end, normal quiet day and 
night models of the lower ionosphere are selected. 

The usual Appleton-Hartree formulas for the elec- 
trical conductivity of a lightly ionized gas of differ- 
ent ionized species are used. Assuming a static mag- 
netic field B to act in the z direction the Cartesian 
form of the tensor permittivity for a harmonic field 
with time dependence exp (—icot) is, 
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The calculations of a t are based on available elec- 
tron and ion density and collision frequency data 
and result for daytime conditions in the curves 
shown in figure 1. The altitude variation of each 
of the three components of the tensor conductivity 
is seen to be nearly exponential through the lower 
part of the ionosphere. 

The propagation in an inhomogeneous, anisotropic 
medium involves, in general, coupled differential 
equations which are difficult to solve. However, 
based on the sharply bounded uniform ionosphere, 
the most pronounced anisotropy effects may be ex- 
pected for propagation along the magnetic equator 
[Wait and Spies, 1960c]. In this particular case 
the surface impedance (required for determination 
of the waveguide mode constants) at the lower edge 
of the anisotropic ionosphere can be derived from 

i The full paper appears in the January 1964 issue of IEEE Transactions on 
Antennas and Propagation. This work was supported by the Office of Naval 
Research under Contract Nonr 3185(00). 
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Figure 1. Components of the daytime tensor conductivity. 



the solution of a differential equation 
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which involves only a single component of the propa- 
gating magnetic field. If the excitation is such that 
no transverse electric field components are present 
(e.g., far fields of a vertical electric dipole of dipole 
moment Idy) the surface impedance in the aniso- 
tropic medium is simply 



Z.=- 
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The differential equation of the required magnetic 
field component can be solved exactly in terms of 
Bessel functions if the relevant components of the 
relative permittivity tensor exhibit an identical ex- 
ponential height variation. A closed form solution 
is obtained for an ionosphere model of two inhomo- 
geneous layers: the lower layer is assumed to be 
isotropic and the upper layer anisotropic. The 
treatment of the isotropic layer is based on solu- 
tions of Galejs [1961a, b, 1962] for an exponential con- 
ductivity-height profile. The assumed exponential 
height variation of the tensor conductivity compo- 
nents does not take into account the gradual leveling- 
off and decrease respectively of these components 
which is noted at heights above 100 km. 

For an arbitrary height variation of the tensor 
conductivity components Z s may be expressed as 



differential equation 
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The new variable V(y) satisfies the Riccati-type 
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Figure 2. Attenuation constants. 
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Equation (6) can be integrated numerically if 
V=V(y u ) is specified at a given altitude y=y u . 
V(y u ) may be determined from the surface im- 
pedance Z s at the upper edge y=y u of the ionosphere 
model, where the differential equation (2) may be 
simplified. Numerical calculations show that y u 
may be selected at y=100 km at daytime and ?/ = 130 
km at nighttime. 

The coupled differential equations are also rela- 
tively simple for propagation transverse to the 
magnetic equator. In this case the coupling becomes 
small, if the vertical variation of the fields is much 
larger than the variation in the direction of propa- 
gation. This condition is obviously satisfied by 
ELF waves. The above differential equations can 
be treated by similar methods. 

The calculated attenuation rates are depicted in 
figure 2. The west-east direction of propagation 
exhibits a lower attenuation constant than the 
east-west direction for/< 1 ,000 c/s. This is contrary 
to the expectations based on a model of a homogeneous 
anisotropic ionosphere. However, this particular 
form of nonreciprocity can be explained with the 
aid of closed form solutions for an assumed exponen- 
tial height variation of the tensor conductivity. 
For /> 1,000 c/s there is a rapid increase in the EW 
attenuation figure, which is principally due to 
absorption in the lower jD-layer, where e x « 1 and 
|e 2 | < C < \l. This increase of attenuation can be 
accounted for after simplifying the differential equa- 
tion (2) for the above values of e x and e 2 and after 
examining its closed form solutions. 

References 

Galejs, J. (Sept. 1961a), Terrestrial extremely-low frequency 
noise spectrum in the presence of exponential ionospheric 
conductivity profiles, J. Geophvs. Res. 66, No. 9, 2787- 
2792. 

Galejs, J. (Nov. 1961b), ELF waves in the presence of ex- 
ponential conductivity profiles, IRE Trans. Ant. Prop. 
AP-9, No. 6, 554-562. 

Galejs, J. (July 1962), A further note on terrestrial extremely- 
low frequency propagation in the presence of an isotropic 
ionosphere with an exponential conductivitv-height pro- 
file, J. Geophvs. Res. 67, No. 7, 2715-2728. 

Wait, J. R. (Feb. 1960a), On the propagation of ELF radio 
waves and the influence of a non-homogeneous ionosphere, 
J. Geophys. Res. 65, 597-600. 

Wait, J. R. (Mar.-Apr. 1960b), Terrestrial propagation of 
VLF radio waves. A theoretical investigation, J. Res. 
NBS 64D (Radio Prop.), No. 2, 153-203. 

Wait, J. R., and K. Spies (Aug. 1960c), Influence of earth 
curvature and the terrestrial magnetic field on VLF prop- 
agation, J. Geophvs. Res. 65, 2325-2331. 

Wait, J. R. (Jan.-Feb. 1962), On the propagation of VLF 
and ELF radio waves when the ionosphere is not sharplv 
grounded. J. Res. NBS 66D (Radio Prop.), No. 1, 53-61. 

(Paper 68D1-322) 



104 



